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ABSTRACT 

The host galaxies of six nearby QSOs are studied on the basis of high resolution HST optical 
images and spatially resolved VLT slit spectra. The gas ionization and velocity are mapped as a 
function of the distance to the central QSO. In the majority of the cases, the QSO significantly 
contributes to the gas ionization in its whole host galaxy, and sometimes even outside. 

Reflection or scattering of the QSO Ha line from remote regions of the galaxy is detected 
in several instances. The line shifts show that, in all cases, the matter responsible for the light 
reflection moves away from the QSO, likely accelerated by its radiation pressure. 

The two faintest QSOs reside in spirals, with some signs of a past gravitational perturbation. 
One of the intermediate luminosity QSOs resides in a massive elliptical containing gas ionized 
(and probably pushed away) by the QSO radiation. The other medium-power object is found in 
a spiral galaxy displaying complex velocity structure, with the central QSO moving with respect 
to the bulge, probably as a result of a galactic collision. The two most powerful objects are 
involved in violent gravitational interactions and one of them has no detected host. 

These results suggest that (1) large-scale phenomena, such as galactic collisions, are closely 
related to the triggering and the feeding of the QSO and (2) once ignited, the QSO has significant 
influence on its large-scale neighborhood (often the whole host and sometimes further away) . 

Subject headings: Quasars - host galaxies - deconvolution - individuals: HE 0306-3301, HE 0354-5500, 
HE 0450-2958, HE 1434-1600, HE 1503+0228, HE 2345-2906 



1. Introduction 

The study of QSOs and their host galaxies has 
become a widespread subject of investigation in 
astrophysics for the past 10 years. Some recent 
developments show that statistical analyses over 
samples of QSOs may lead to fruitful conclusions 
about the morphology, the gaseous and stellar con- 
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tent of the host galaxies, and about the mass of 
the central black hole, the QSO ignition and fu- 
eling, etc. However, no clear scheme has emerged 
yet to account for the diversity of QSO and host 
galaxies properties and the debate on their inter- 
relations is still widely open. 

Previous studies of QSO host galaxie s were car- 
ried out using either spectros c opic dat a ( Hughe s et al 



2000; Baldwin et al. 1981; Letawe et all 



Vanden Berk et al.. I 120061 ; iKewlev et al 



or high reso l ution images (Bahca ll et al 



2007 



2006) 



1997 



Dunlop et all 120031; iFlovd et alj |2004) . None of 



them used both spectroscopy and high resolution 
imaging. 

This kind of study was carried out for the 



bright QSO HE 0450-2958 (|Magain et all 12005') 



showing a strongly distorted galaxy, about 7 kpc 
away from the QSO. Previous studies (jBovce et al 
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1996; ICanalizo et aL 2001 ). assumed that it prob- 
ably experienced a collision with the QSO host 
galaxy. A careful processing of the HST images 
allowed to separate the point source from the ex- 
tended objects, with the surprising result that no 
host galaxy could be detected around the QSO. 
Only a compact 'blob' next to the QSO could be 
found, whose spectrum was shown to consist in a 
cloud of gas ionized by the QSO radiation, with no 
trace of a stellar component. The nature of HE 
0450-2958 is still controversial: is there no host 
galaxy at all, or could it be faint enough to escape 
detection? Which formation mechanism could ac- 
count for such a special configuration? Whatever 
the answer, the study of HE 0450-2958 is a beau- 
tiful example of the nice complementarity between 
spectroscopy and imaging. 

The aim of the present paper is to extend this 
kind of analysis to the six Type 1 QSOs for which 
we have both high resolution images obtained with 
the Hubble Space Telescope (HST) and slit spectra 
using F ORS1 on the ESQ Very Large Telescope 
(VLT) (|Letawe et al.ll2007t ). The six QSOs are: 
HE 0306-3301, HE 0354-5500, HE 0450-2958, HE 
1434-1600, HE 1503+0228 and HE 2345-2906. 

Section 2 gives an overview of the observational 
data and of the main characteristics of the reduc- 
tion process. In Sect. 3, the post-processing of the 
HST images is explained, while Sect. 4 describes 
the extraction of the spatially resolved host spec- 
tra. In Sect. 5, the methods used to analyze si- 
multaneously the images and the spectra and to 
infer the host properties are explained. Section 
6 presents the results on individual objects, while 
Sect. 7 provides a general discussion and some con- 
clusions on the QSO host properties. Through- 
out the paper, we adopt the following cosmology: 
H = 65 km s -1 Mpc-\ O m = 0.3 and fl x = 0.7. 

2. Observations and reduction 
2.1. The sample 

The six QSOs form a subsample of the 20 bright 
low redshift QSOs (z < 0.35, M v < -23) from the 
Hamb urg-ESO Survey analysed by iLetawe et al 



(|2007j ) for which we also have high resolution im- 
ages obtained with HST. The subsample was se- 
lected on the basis of their spectra with preference 
to those revealing large amount of gas and interest- 
ing features, which introduces unavoidable biases. 



Moreover, It is obviously too small to attempt any 
statistical analysis. Thus, the present analysis can 
be seen as (1) a demonstration of the interest of 
combining high resolution images and spatially re- 
solved spectra to study the QSO hosts and (2) the 
focus on interesting cases, which can shed light on 
the relations between the QSOs and their hosts 
galaxies or, more generally, their environment. 

2.2. Images 

The images have been obtained in September- 
October 2004 and April 2005 with the Advanced 
Camera for Surveys/High Resolution Channel on- 
board HST with the broad V-band filter F606W. 
The ACS/HRC detector consists of a 1024 x 1024 
pixels CCD, with - 0.028" x 0.025" per pixel, cov- 
ering a 29" x 25" field of view. 

For each QSO, six exposures, taken during one 
HST orbit, are available: three short unsaturated 
ones, in order to be able to study the host prop- 
erties very close to the central QSO, and three 
long ones, saturated around the center, but pro- 
viding a good S/N in the lower surface brightness 
regions. Small offsets are applied to ease cosmic 
rays and bad pixels removal (5 pixels in each di- 
rection between successive exposures). Moreover, 
a PSF star is observed after each QSO (Fig. [T]). 
Table [1] summarizes the main QSO characteristics 
and the characteristics of the HST observations. 

The data are provided under different forms, 
along with auxiliary files: 

- the observation itself (available in raw, flat- 
fielded, or drizzled form) 

- a flag image indicating hot, warm or bad pixels, 
and cosmic rays hits 

- a standard deviation image. 

The drizzled data are the individual raw images 
corrected from optical distorsions inherent to the 
ACS instrument and combined into a single im- 
age. This correction uses interpolations. We avoid 
their use in order to lower the number of manipula- 
tions, in particular interpolations which introduce 
correlated noise in adjacent pixels. We rather cor- 
rect for distorsion at the very end of the process, 
on the subsampled deconvolved images. Our data 
processing method (MCS, see below) allows both 
to simultaneously deconvolve several dithered im- 
ages of the same field and to oversample the data 
by a factor of 2. It thus performs its own optimal 
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Table 1: Observational char acteristics of the sa mple (redshift, absolute V-band magnitude of the QSO My 
from the VLT observations (jLetawe et al.l 120071 ). short and long exposure times and position angle (PA) of 
the observation, counted positive East of North. 

Name Redshift My (QSO) Short E.T. Long E.T. PA 









(sec) 


(sec) 


(°) 


HE2345-2906 


0.223 


-22.97 


40 


340 


177.8 


PSF 






10 


30 




HE1503+0228 


0.135 


-23.43 


40 


330 


36.8 


PSF 






10 


40 




HE1434-1600 


0.144 


-23.82 


30 


330 


288 


PSF 






5 


40 




HE0306-3301 


0.247 


-24.81 


30 


330 


110 


PSF 






15 


45 




HE0354-5500 


0.267 


-25.45 


40 


400 


258 


PSF 






5 


40 




HE0450-2958 


0.285 


-25.82 


30 


330 


105.6 


PSF 






G 


45 






drizzling. The image reduction itself is performed 
on the flatfielded data using the Pyraf software. A 
standard deviation map is created using, for every 
pixel i of an observation, the formula 



(1) 



where Aj is the signal in pixel i, g the CCD gain, 
and r the background (readout) noise. Then we 
create an inverse sigma image. This allows to give 
stricly nil weight to bad or saturated pixels and 
cosmic ray hits, simply by attributing them a zero 
value on the inverse sigma image. This corrected 
inverse "noise-image" is used for the MCS process- 
ing. 

2.3. Spectra 

The spectroscopic data were obtained with 
the ESO VLT/UT1 in April and November 2000 
with FORS1 in Multi Object Spectroscopy (MOS) 
mode, which consists in a set of l"-slits with a 
fixed length of 19". Every QSO was observed 
with 3 different filters, covering a total restframe 
wavelength band of 5000A, from 3500A to 8500A. 
One slit was centered on the QSO and the other 
ones on PSF stars and surrounding objects. 

The reduction and deconvolution meth ods are 
presented in detail in lLetawe et alJ (|2007l ). result- 



galaxy. This spatial separation of the point-like 
component from the diffuse one relies on the hy- 
pothesis that the latter contains no significant 
structure narrower than the fixed and finite res- 
olution of the deconvolved spectrum. In order to 
regularize the solution of this inverse problem, one 
has to apply a smoothing constraint on the diffuse 
component. 

Such inevitable smoothing decreases the spatial 
resolution of the host galaxy spectrum, especially 
in cases where the host S/N is low (the higher 
the S/N, the better the spatial resolution of the 
host spectrum). To avoid this variable-resolution 
proble m, we proceed i n a sli ghtly different manner 
than in lLetawe et al.l (|2007t) . Instead of relying on 
the host spectrum as given by the deconvolution 
process, we start from the original (reduced) spec- 
tra and subtract the two-dimensional QSO-alone 
spectrum given by deconvolution. Thus, we get a 
2D host galaxy spectrum with a spatial resolution 
constant along the slit and equal to that of the 
original data. 

3. Image processing 

The HST/ACS images are pr ocessed with 



the M CS deconvolution algorithm (M again et al 



ing in spatially deconvolved spectra with the QSO 
spectrum separated from the spectrum of its host 



Il998h . which is particularly well suited for separat- 
ing point sources from diffuse ones and, as such, 
very well adapted to the study of QSO host galax- 
ies. Indeed, as demonstrated earlier in the case of 
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HE 0450-2958 (jMagain et alj|2005l) . the extended 
wings of the HST PSF may cause the light from 
the bright QSO to significantly contaminate the 
image of the nearby souces. Important compo- 
nents can be hidden, such as "the blob" in the 
aforementioned case. 

The principle of the MCS method is to produce 
images with a resolution and sampling fixed by the 
user, provided that they conform to the sampling 
theorem. In our case, the deconvolved image has a 
pixel size two times smaller than the original HST 
images and a PSF of Gaussian shape, with a Full 
Width at Half Maximum (FWHM) of two final 
pixels. The point sources are explicitely separated 
from the diffuse components, which allows to get 
an image of the host galaxy uncontaminated by 
the QSO (see Fig. |2J|. 

3.1. PSF construction 

A crucial point in the deconvolution process is 
the PSF construction: the more accurate the PSF, 
the better the deconvolution. In order to con- 
struct the be st possible PS F, we use the method 
explained in iLetawd (|2006f ). Here are its main 
steps: 

(1) With the Tiny Tim softwanQ, we create a 
synthetic PSF corresponding to the actual observ- 
ing conditions. Let us call it TT. We then com- 
pute the kernel which allows to transform the TT 
image into a gaussian r of 2 pixels FWHM. Let us 
call this kernel Tiny Dec. It obeys the equation: 

TT(x) = TinyDec(x) * r(x), 

where a; is a 2D vector representing the posi- 
tion in the image and * stands for the convolution 
operator. 

(2) TT is a model of the actual PSF. However, 
it is not accurate enough to provide a trustwor- 
thy separation of the QSO from its host galaxy. 
We improve it by comparing TinyDec(x) *r(x) to 
the observation t(x) of a point source (one of the 
reduced PSF star observations). In order to im- 
prove the fit, we add a numerical component F{x) 
to the model TinyDec(x) until the observed PSF 
star image is reproduced with satisfying accuracy: 

t(x) = (TinyDec(x) + F(x)) * r(x). (2) 



Among all the PSF stars, we choose the one which 
is best suited for a given QSO observation. As 
the HST PSF depends on a number of param- 
eters (such as observed target spectral distribu- 
tion, image position in the focal plane, position 
of the satellite on its earth orbit,. . . ) we test a 
variety of QSO/PSF combinations and choose the 
one which gives the smallest residuals in the point 
source modelling. 

3.2. Simultaneous deconvolution 

Once the PSF is chosen, we proceed to a si- 
multaneous deconvolution of the six individual im- 
ages, in order to reach the best possible resolution 
and S/N ratio. As mentioned above, the MCS 
method is able to separate a pointlike source from 
a diffuse background, i.e., in our context, the QSO 
from its host galaxy. 

We find that the best results are obtained in 
the following way: 

(1) The position and intensity of the point source 
in each individual image are determined through 
deconvolution by ensuring that the image with the 
QSO point source subtracted contains a minimum 
amount of structure from the point source PSF 
(Airy rings, spikes,. . . ). In this way, we obtain in- 
dividual host galaxy images. 

(2) In a second phase, the MCS program is run 
on the six host galaxy images, which are simulta- 
neously deconvolved, and a model is found which 
best represents this full set of observations. The 
model is separated into two components: an an- 
alytical model accounting for the galaxy bulge, 
and a numerical background image accounting for 
all supplementary features (spiral arms, neigh- 
bors,. . . ). 

The analytical model is chosen to follow a "de 
Vaucouleurs" law with elliptical isophotes, i.e. an 
intensity profile: 



I(x) = I(xi,x 2 ) = he 



{ Ax±+Bx 2 + Cx 



1X 2 ) 



(3) 



J -Krist J., 

http: / / www.stsci.edu / software / tinytim / tinytim . html 



where x\ and X2 are the two coordinates relative 
to the bulge center. The free parameters Io, A, B 
and C are determined by the program. 

An example of the whole process is illustrated 
in Fig. 1 for the QSO HE0354-5500. One sees 
that the residuals are very good, except around the 
very center. The reason for this misfit in the cen- 
tral region is that, in long exposures, both the PSF 
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star and the QSO are saturated. As the weights of 
the pixels in the saturated regions have been set 
to zero, they do not constrain the fit. 

4. Spatially resolved spectra 

By superimposing the projection of the VLT 
slit onto the HST host galaxy images, we can iden- 
tify which parts of the galaxy contribute to each 
part of the spatially resolved spectra. Of course, 
we have to take into account the seeing condi- 
tions during the ground-based spectroscopic ob- 
servations. This is done by convolving the host 
galaxy images with a 2D Moffat function derived 
from the ID profile of the PSF stars observed si- 
multaneously with the QSOs, and assuming circu- 
lar symmetry of the PSF, which is accurate enough 
for the present purpose (see Fig. [3]). Even with the 
degraded resolution, it is still meaningful to sepa- 
rate the spectra into spatially resolved elements of 
different parts of the galaxy, such as spiral arms, 
bulge or halo. Several examples of extraction are 
given in Fig. [4] for the H/3-[OIII] spectral region. 
These spectra are analyzed in Sect|6] 

5. Spectral analysis 

Before any analysis, we correct the spatially re- 
solved spectra for reddening, assuming a theoret- 
ical Balmer decrement Ha/H/3= 3.1 if the region 
is ionized by the nucleus, and Ha/H/3= 2.85 if it 
is ioniz ed by stars, using the W hitford reddening 
curves dMiller fc Mathewslll972f) . The forthcom- 
ing detailed analysis focuses on the gaseous prop- 
erties of the galaxies only, as the S/N for absorp- 
tion lines is too weak. 

5.1. Diagnostic diagrams 

A diagnostic diagram uses line ratios to distin- 
guish between the main ionization sources of the 
emitting gas, i.e. star formation, ionization by 
the QSO radiation o r by shocks, as suggested by 



Baldwin ct al. | (|l98lf) and lVeilleux fc Osterbrock I 
1987 ). These diagrams were revised bv lKewlev et all 
2006), who added a composite region, that 



is, a region in the diagram where the ioniza- 
tion process is due to a mixture star forma- 
tion and QSO radiation. The fluxes in the 
main emission features are measured from the 
spectra of the different galactic regions, as de- 
fined in the previous section. We then compute 



the emission line ratios [OIII](5007)/H/3(4861), 
[NII](6583)/Ha(6563) and [Oil] (3727) /[OIII] (5007) 
for each region and plot them in diagnostic dia- 
grams. 

Each region is differently affected by uncertainties 
in the measurements and so should have its own 
error bar. For clarity, when these uncertainties 
do not vary much from region to region, we only 
put one global error bar, which is representative 
of the different regions. Otherwise, specific error 
bars are displayed for each region. 
Up to now, such diagnostic diagrams had only 
been used for Type 2 AGNs (i.e. narrow-line 
AGNs) and for the central regions of Type 1 
(broad and narrow-line) AGNs. As far as we 
know, this is the first time that the hosts of Type 1 
QSOs are analyzed with such diagnostic diagrams, 
except in iLetawe et al. ( 2007 ) where the diagnos- 
tic diagrams are constructed from the integrated 
spectra. 

5.2. Radial velocity curves 

A very useful tool in the investigation of a QSO 
host is the radial velocity curve. We construct 
such a curve by comparing the wavelength shifts 
of some emission lines (e.g., Ha or [OIII]) for every 
pixel along the slit. Then we transform the shifts 
in relative velocities with respect to the central 
QSO. An analysis of the velocity distribution of 
the host relative to the QSO is often fruitful to 
understand global tendancies in the dynamics and 
to model the mass distribution in the host. All 
the radial ve locity curves and t heir analysis can 



be found in (jLetawe et al.l 120071 ) . Some of them 



are discussed in Sect|7J] 

5.3. Line fits 

The Ha(6563A) line is partially blended with 
the ncig hboring [Nil] lines (at A6548A and A6583A). 
In order to get accurate values for the line posi- 
tions, widths and intensities, we model them as a 
sum of gaussian profiles. 

The two [Nil] lines have a theoretical intensity 
ratio of 

I([NII] (6583)) 



3. 



(4) 



I([NII] (6548)) 

according to lOsterbrockl (1989). This ratio is kept 
fixed in the fit as it only depends on the atomic 
structure and not on the astrophysical environ- 
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ment. Moreover, both widths are set equal be- 
cause the [Nil] emissions arise from the same re- 
gion, characterized by the same temperature and 
velocity distribution. 

These emission lines are generally superim- 
posed on a stellar spectrum, consisting in a con- 
tinuum plus absorption lines. In all the cases an- 
alyzed, the emission lines are so strong compared 
to the continuum that the latter can be locally 
approximated by a straight line, thus neglecting 
stellar absorptions and continuum curvature. 

Taking these constraints and assumptions into 
account, we fit the spectrum around Ha using a 
sum of a linear continuum and three gaussians 
with two restrictions on the parameters: the in- 
tensity ratio between the [Nil] lines as well as their 
widths. Some of the spectra are very well fitted 
by this method (see Fig. [5] for an illustration), 
but this fit turns out to be unsuccessful for some 
regions located far away from the nucleus in HE 
0354-5500, HE 0450-2958, HE 1503+0228 and HE 
0306-3301. 

The three-gaussians model described above is 
unable to fit the data around the Ha line in those 
regions, and an additional broader component is 
necessary. Adding a fourth gaussian generally al- 
lows to get an acceptable fit. 

This fourth gaussian is much broader than the 
other components and its width is similar to that 
of Ha in the QSO spectrum. We further test this 
similarity by assuming that the broad component 
has exactly the same profile as the QSO Ha line. 
Under such an assumption, we proceed to fit the 
spectrum in two steps: 

• Step 1: we fit the Ha line of the central QSO 
spectrum by a sum of two gaussians (see Fig. 

• Step 2: we fit the spectra of the host regions 
by a sum of 5 gaussians: 3 gaussians to ac- 
count for the narrow Ha and [Nil] lines with 
the constraints described above. The two re- 
maining gaussians are obtained from the fit 
of the QSO Ha line in the Step 1. Their pa- 
rameters are constrained to have the same 
shape as the QSO Ha line: the intensity ra- 
tio, the shift between the centers and the 
widths are fixed. So, the only free parame- 
ters are the global intensity and the center 
the line. 



However, we have to make sure that this broad 
line is not due to an imperfect subtraction of the 
QSO spectrum and that it indeed corresponds to 
a real feature. 

The major source of error in deconvolution is 
the quality of the PSF. However, we have tested 
that changing the PSF within reasonable limits 
does not remove this broad component. More- 
over, if the broad line did arise from an imperfect 
subtraction of the QSO spectrum, it would most 
likely appear in all the regions of the host, because 
each wavelength is deconvolved independently in 
the spatial direction. But this broad component 
only appears in some regions, always located quite 
far from the nucleus. 

In most instances, the broad component ap- 
pearing in the host spectrum is redshiftcd in wave- 
length with respect to the QSO emission line, 
while an improper subtraction of the QSO spec- 
trum would have most likely left a line at the same 
wavelength as in the central QSO. 

We have also tested if the use of Lorentzian pro- 
files could improve the fit to the data, as one might 
expect that the broader Lorentzian wings would 
weaken the need for an additional broad Ha com- 
ponent. However, it turns out that a combina- 
tion of Lorentzians does not provide an adequate 
fit. Lorentzian profiles can be used to fit the spec- 
tra of some AG N when they have broader wings 
(|Veilleuxj|l991l ). but in our case, the spectra are 
coming from different parts of the host galaxy and 
not from the AGN itself (the QSO spectrum has 
been removed) . In such a case, the emissions lines 
are expected to be gaussian-shaped and, indeed, 
gaussians are found to provide satisfactory fits in 
the majority of the cases (with the notable excep- 
tion of those - remote - regions presenting a broad 
component). 

The improvement of the residuals of the fit 
when we add the QSO-like component in the case 
of HE 0354-5500, is shown in Fig. 

5.4. About the origin of the broad Ha line 

Which mechanism causes broad lines to be 
present in the spectrum of the host galaxy? The 
same kind of phen omenon is known in the radio 
galaxy Cygnus A ( Ogle et al~l 1997), the Se yfert 



2 galaxies N GC 1068 (ISimpson et l2002h and 



NGC 4388 (IShields fc Filippenko I Il996l) . and in 
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luminous infrared galaxies ( Tran et al. 2000). In 
all those cases, the nucleus light is obscured or un- 
observed but emission coming from the nucleus is 
either reflected by dust or scattered by free elec- 
trons (for simplicity, we shall call it the reflected 
light). This reflection (or scattering) causes the 
light to be polarized. Spectro-polarimetry is then 
a good way to observe reflection of the nuclear 
emission in unresolved objects, and the analysis 
of the polarized reflected light allows to study in- 
directly some properties of the nucleus, such as 
its orientation or luminosity, while not being able 
to observe it dire ctly because of dust obscuration 
(|Antonuccilll993h . 



in a resolved region of the host by the for- 
mula 



Av 



,(1 — cos 



A few cases of broad Mgll or H/3 line reflec- 
tion, observed with direct spatially resolved spec- 
tros copy, are reported in Dev fc Spinrad ( 19961 ) 
and ICid Fernandes et alj (|2004l ). but, once more, 
only for Type 2 objects (for which the contrast 
host/QSO is much favourable than in Type 1 
QSOs). In our case, only the Ha line is strong 
enough to be detected by reflection in some re- 
gions of the QSO environment. 

In our sample of Type 1 QSOs, the similar- 
ity between the observed broad Ha lines in the 
host galaxy and the corresponding QSO Ha lines 
clearly suggests that reflection or scattering might 
be the explanation. The main characteristics of 
the reflected broad lines are listed in Table [2] which 
gives: 

1. The relative velocity shift of the reflected 
line with respect to the QSO emission line. 
A shift between the centers of the reflected 
Ha line and the QSO emission line clearly 
appears in some regions. In Tabled all those 
shifts have been transformed in relative ve- 
locity Av. Error bars on the measured shifts 
are estimated by shifting the center of the 
broad Ha until the fit becomes unaccept- 
able. In a basic case of reflection, a shift 
between the emitted and the reflected light 
is interpreted in terms of motion of the re- 
flecting medium. Thus, a shifted broad line 
reveals the presence of some material in mo- 
tion with respect to the QSO. The fact that 
all shifts are positive shows that the reflect- 
ing medium is moving away from the QSO. 
The observed shift Av is linked to the expan- 
sion velocity v oxp of the scattering particles 



where 9 is the angle between the direction 
of motion of the reflecting medium and the 
line of sight. 

Consequently, the measured shift must be 
seen as a lower limit to the expansion veloc- 
ity of the reflecting medium with respect to 
the QSO. 

2. The flux ratio R = , ^f lected ^ between the 

Fqso 

flux of the reflected Ha broad line and the 
flux of the Ha line in the QSO spectrum. 

3. The equivalent width (EW) of the reflected 
line compared to the QSO emission line. If 
we assume that the medium reflects both 
the Ha line and the underlying continuum, 
and that a stellar continuum may also be 
present, the reflecting regions must have an 
EW equal (if there is no stellar continuum) 
or smaller (in absolute value) than the QSO 
emission line. This is indeed the case. 

4. The ratio § a , evaluated on the galactic nar- 
row component of Ha for the region contain- 
ing the broad line, reveals the presence of 
dust, as reddening increases this ratio above 
the nominal value of 3.1 or 2.85, depending 
on the excitation mechanism. 

5. The deduced electron temperature and an 
upper limit on this temperature. If the scat- 
tering is due to electrons, the reflected light 
must be broadened according to the random 
motion of the scattering electrons. In prac- 
tice, we convolve the QSO Ha line by a gaus- 
sian representing the electron velocity distri- 
bution (corresponding to a given tempera- 
ture) and use this new profile instead of the 
QSO Ha line in the 5-gaussians fit. The best 
fit gives an estimate of the electron temper- 
ature T e - . Increasing the electron tempera- 
ture until the fit becomes unacceptable sets 
an upper limit T max . 

All regions containing a reflected line are faint 
in the HST images (except in HE 1503+0228), re- 
vealing a weak (if any) stellar continuum. Indeed, 
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Table 2: Characteristics of the reflected brad Ha line in the spectrum of the host galaxies. 



Object 



Region Rel. vel. 
(km s' 1 ) 



Flux ratio R 



HE 1503+0228 
HE 1503+0228 
HE 0306-3301 
HE 0354-5500 
HE 0450-2958 



A 
E 
F' 
E 
B 



effected Hq/QSO He 
\EW\ 



c (10 5 if) 



+620t^° 
+458 ± 180 
+745+ 250 

^ — 280 



3.09 ± 0.09 10~ J 
3.36 ± 0.10 10~ 3 
1.03 ± 0.05 10~ 3 
4.06 ± 0.12 10~ 3 
9.33 ± 0.20 10~ 3 



71/345 
98/345 
280/365 
117/116 
422/558 



5.38 
4.28 



3.66 
3.49 



2.27 - 6 
5.55 - 12 + 
0.73 - 1.86 
1.02 - 4 
- 1.35 



one might imagine thet the reflected Ha line may 
be present in many regions, but only observable 
when the stellar continuum is weak enough, that 
is, in low luminosity regions of the images. 

It is hard to determine whether the phe- 
nomenon at work is mostly reflection by dust or 
scattering by free electrons. First, the |p ra- 
tios in Table [2] indicate that all regions showing 
a broad Ha line contain some dust. Second, all 
fits with zero broadening can be considered ac- 
ceptable, which shows that dust reflection is a 
viable hypothesis. On the other hand, in all ob- 
jects except HE0450-2958, the best fit is obtained 
with a line having a non-zero broadening, which 
points towards electron scattering. Moreover, the 
temperatures obtained under the assumption of 
electron scattering appear reasonable and consis- 
tent with the estimates of Dev fc Spinrad ( 1996f ) 
(T e - w 3.10 5 .fO. Finally, all the spectra of these 
regions show high ionization lines, which indicates 
that free electrons are present. It is thus very 
difficult to decide between the two different alter- 
natives. Both may contribute to various degrees 
in different objects. 

6. Results on individual objects 

The main characteristics of the host galaxies^, 
as derived from the HST images, are summarized 
in Table [3l The individual cases are discussed be- 
low. We start by the least powerful nuclei, HE 
1503+0228 and HE 2345-2906, and proceed in or- 
der of increasing nuclear luminosity. 

6.1. HE 1503+0228 

The a nalysis of this object h as already been pre- 
sented in lCourbin et alj (|2002h . The HST images 



confirm their conclusions, but bring significant ad- 
ditional information. 

At first sight, the host of HE 1503+0228 seems 
to be a quite mundane spiral galaxy (Fig. [8]). 
However, the deconvolved image reveals a faint, 
extended component with approximately circular 
symmetry, surrounding the whole system at ss 6.5 
kpc from the center. 

As can be seen from Fig. [8j the different regions 
are clustered in the same area of the diagnostic 
diagram, indicating a mix of ionization by shocks 
and by stars. However, this faint extended com- 
ponent (regions A and E) differs from the bulge 
and from the inner arms in that it shows a broad 
Ha component superimposed on the narrow emis- 
sion lines and on a significant continuum. Among 
all the broad reflected Halines, those detected in 
HE 1503+0228 have the smallest EWs. This sug- 
gests that the contribution of the stellar contin- 
uum to the total spectrum is significant. This 
conclusion is further supported by the fact that 
the spectrum around Ha cannot be properly fit- 
ted with a 5-gaussian model. Taking into account 
a stellar Ha in absorption allows to improve the 
fit (see Fig. [5]). The x 2 minimization leads to 
an EW of the absorption line of 11.5 A for re- 
gion A and 16.5 A for region E, which is larger 
than gener ally found in inte g rated g alactic spectra 



(see, e.g., iKennicutt et al.l ( 19921) . Tresse et al 
( 19961 )). However, the S/N of our spectrum is 



The bulge ellipticity is defined in terms of our analytical 
model by £ = , D , , where 8 = 0,5 arctan( „ c . ). 



rather low in regions A and E, so that the uncer- 
tainties are quite large (the one sigma uncertainty 
on the intensity of the absorption is about 35%). 
The large EW of the absorption line suggests that 
the stellar spectrum in these external regions is 
dominated by a young stellar population. This is 
consistent with the conclusion drawn from the di- 
agnostic diagram, which shows a mix of ionization 
by shocks and by stellar light, and with the large 
Balmer decrement measured in these regions (Ta- 
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Table 3: Global properties of the host galaxies. Morphology, nucleus to host luminosity ratio (N/H), bulge 
ellipticity (bulge ell.) and orientation of the main axis counted positive East of North (9 in degrees). 
HE0450-2958 is not listed, as no host was detected. 



Name 




Morphology 


N/H 


bulge ell. 


9 


HE 1503- 


f0228 


Spiral 


1.9361 


0.0116 


10 


HE 2345- 


-2906 


Spiral 


0.9039 


0.0294 


-43 


HE 1434- 


-1600 


Ellipt. 


1.7833 


0.1698 


-29 


HE 0306- 


-3301 


Spiral 


7.4079 


0.3921 


37 


HE 0354- 


-5500 


Spiral 


2.1568 


0.2534 


-101 



ble [2]) . Indeed, if significant stellar Ha is present 
in absorption, it should also be seen in H/?. Such a 
stellar absorption should introduce a larger under- 
estimate of the relatively weak H/3 narrow emis- 
sion line, when compared to the much stronger Ha 
emission (thus explaining the high Balmer decre- 
ment measured in these regions). 

The presence of a broad reflected line in the 
outer regions A-E and not in the inner ones B- 
D raises the following problem. Since these re- 
gions show otherwise similar spectra, with about 
the same ionization level, one could wonder how 
the QSO light can cross the inner spiral arms with- 
out being absorbed or reflected, to reach the outer 
regions, and be reflected there. The most natural 
explanation seems to be that these outer reflect- 
ing regions are not located in the galactic plane. 
This is supported by the deconvolved HST image, 
which shows that, while the inner region, contain- 
ing the spiral arms, has a significant ellipticity, 
the outer part is much more circular. If it was lo- 
cated in the same plane, and had a circular sym- 
metry, the inclination angle of the galactic plane 
should introduce the same apparent ellipticity as 
for the inner region. Its circular morphology sug- 
gests that it is either in a different plane, close to 
the plane of the sky, or that it has a roughly spher- 
ical symmetry (e.g. a halo populated with young 
stars). These peculiarities are supported by the 
radial velocity curve (Fig[5]), which shows a rather 
sharp drop in the outer E region (the spectrum is 
more noisy in the A region). This drop could be 
interpreted as due either to a rotation in a plane 
closer to the plane of the sky or to a geometry 
closer to spherical. 

In the case of non active galaxies, outer rings sur- 
rounding spirals are supposed to be created either 
by gravitational interactions (merger, collision, or 
even tidal interactions by a neighboring galaxy), 



or by resonance effects with a bar, if any (jButa 



199(1 iButa et all Il998l Il999h . In our case, the 



peculiar geometry and spectra of HE 1503+0228's 
outer ring suggest that some gravitational interac- 
tion probably happened in the past (e.g. a minor 
merger or a tidal interaction with a passer-by). 
This past interaction might also be related to the 
ignition of the QSO activity. 

6.2. HE 2345 2906 

As seen from Fig. [TU1 the host is a barred spi- 
ral, with a prominent 8 kpc-long bar. Two spi- 
ral arms originate from the extremities of the bar 
and surround the close-to-spherical bulge. A wider 
and fainter arm, located 8 kpc from the center, 
surrounds the whole system. The diagnostic di- 
agrams show that the emission regions lie in the 
composite region of the diagram, which means that 
the gas ionization is probably due to a mix of 
shocks, stellar light and QSO radiation. 

A bar in spiral host galaxies is often consid- 
ered as an efficient means of funnelling matter 
into the central r e gions and feeding the AGN (e.g. 
Sellwood I (|2000h : ICrenshawl (|2003l )V As shown 
in Fig. no broad reflected component is de- 
tected. Moreover, all the spatially resolved spectra 
are typical of a classical spiral galaxy. 

6.3. HE 1434 1600 

The spect r a of t his object were analyzed in 



Letawe et al.l (|2004f ) . and their results are con- 
firmed by the present study. 

The deconvolved image (Fig. [TTj) shows that, 
apart from the central QSO, the system is made 
of at least three components: (1) an elliptical host 
completely taken into account by our analytical 
model (whose central part corresponds to region 
B), (2) a neighboring elliptical galaxy, and (3) 
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some filamentary structures (regions A, C and D) 
on both sides of the center and ~ 3 to 8 kpc far 
away from it. The ionization diagrams of Fig. 1111 
show that the filaments consist of gas highly ion- 
ized by the QSO. The central region B lies closer 
to the transition between ionization by shocks and 
by the AGN. 

The analysis of t h e rad ial velocity curve (Fig. 
15 of Letawe et al. ( 20041 )). enlightened by the 
ACS image, shows that the filaments are moving 
with respect to the central QSO, at a speed of 
at least 150 km/s (which is the measured radial 
velocity), and slightly increasing outwards. They 
most probably consist of gas blown away and ac- 
celerated by the QSO radiation. 

The presence of extended regions of highly ion- 
ized gas in an elliptical galaxy is unusual. It 
might have been brought there by a past inter- 
action (maybe related to the QSO ignition), and 
its presence is revealed to us only thanks to the 
QSO strong ionization field. The bipolar struc- 
ture of the gas is probably related to an orienta- 
tion effect of the QSO. Radio o bservations of HE 
1434-1600 (jCondon et alJ Il998h reveal two lobes 
extending up to 980kpc, North and South of the 
QSO, which is roughly the same orientation as the 
gaseous filaments (Fig. [TTj) . Thus, the ionized gas 
might only be a part of the whole gas present the 
host galaxy. 

6.4. HE 0306 3301 

Figure [12] shows the deconvolved image of the 
host galaxy of HE 0306 — 3301. The deconvolution 
of this particular case has been harder to handle 
as the QSO is highly saturated and the central 
region affected by the loss of linearity extends to a 
larger area. This explains the residual spikes and 
apparent gaps around the very center. 

At first glance, the host looks like a rather mun- 
dane spiral galaxy. However, deeper scrutiny re- 
veals that this first impression is far from telling 
the whole story. 

First, the spiral arms look asymmetric and are 
more prominent in the N-W part of the host. In 
this region, a compact source (that we shall call 
'the spot') is found in (or close to) one of the spiral 
arms, at about 5 kpc from the QSO. A close-up 
on that region is shown in Fig. [121 It is well fitted 
by a point source, plus a faint extension in the 



direction of the QSO. In the V-band, the spot is 
about 1000 times fainter than the central QSO. 

The diagnostic diagrams presented in Fig. [12] 
show that all regions contain gas whose ionization 
level is in the transition region between stellar and 
AGN (or shocks). There is a tendency for the 
ionization to increase from the S-E to the N-W, 
i.e. when crossing the galaxy in the direction of 
the spot and further away. 

The VLT spectra are not deep enough to al- 
low extracting the spectrum of the spot. How- 
ever, they reveal a very unusual radial velocity 
curve over the whole galaxy (Fig. [T2|) . The N- 
E side of the host has a rather constant velocity of 
~ 250km s _1 relative to the central QSO. How- 
ever, when approaching the QSO position, in the 
inner 1 arcsec, the emission lines become resolved 
into two components. We have fitted two gaus- 
sians to separate these two components and find 
that, in all three main emission lines, i.e. Ha, 
[Nil] and [OIII], the highest velocity component 
remains at about the same velocity as in the outer 
regions. Its position in the diagnostic diagrams 
shows a mix of ionization sources, as found in the 
other regions of the galaxy. The lowest velocity 
component displays radial velocities approaching 
zero in the QSO rest frame, as well as an ionization 
degree typical of star-forming regions. 

The depth of the spectra does not allow to build 
a radial velocity curve based on the absorption 
lines. However, an average absorption lines radial 
velocity is determined on each side of the QSO by 
correlating the host spectrum with a template spi- 
ral galaxy spectrum in the spectral region between 
the [Oil] and H/3 lines, which is devoid of measur- 
able emission lines. We obtain +175 km s -1 for 
the S-E side and +80 km s _1 on the N-W side, to- 
wards the spot. The error bars are of the order of 
50 km s _1 . This means that the average velocity 
of the bulge stars is shifted by ~ 130 km s _1 with 
respect to the QSO. 

Velocity shifts between different QSO emis- 
sion lines are common, QSO lines being generally 
blueshifted with respect to their host galaxy rest- 
frame. One can thus wonder how to obtain a reli- 
able measurement of the QSO redshift . This ques- 
tion was investigated bv lLetawe et al.l ( 20071 ). who 
compared the velocity shifts of different QSO emis- 
sion lines with the host galaxy radial velocity, in 
five cases of spiral hosts with regular morphology 
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and symmetrical rotation curves. In such cases, 
the host redshift can be determined accurately 
and without ambiguity. They found that, in these 
symmetrical cases, the velocity measured from the 
tip of the QSO Ha and H/3 emission lines agrees 
with the host velocity, while the bulk of these lines 
indicates a blueshift (related to the well-known ex- 
tended blue wing). Narrow QSO emission lines are 
also generally blueshifted. 



Following iLetawe et al.l (|2007l ). we have used 
the tip of the Ha and H/3 emission lines to es- 
timate th e QSO redshi f t. Wh en compared to the 
results of ILetawe et al.l ( 2007 ). which show an av- 
erage velocity shift compatible with zero and a 
scatter of about 20 km s _1 , the 130 km s _1 differ- 
ence found in this case is clearly an outlier. The 
velocity difference we measure between the QSO 
and the bulge of its host galaxy can be considered 
secure at the 6<r level. 

Considering now the radial velocity curve on 
the N-W side of the QSO, Fig. [T2l shows that all 
three emission lines indicate roughly the same ve- 
locity in the inner 1 arcsec, ~ —50 km s _1 with re- 
spect to the QSO, or ~ —180km s _1 with respect 
to the average velocity of the stars. Then, the 
radial velocities of the three emission lines start 
to increase at different rates. The [OIII] emission 
line increases first, followed by [Nil] and then by 
Ha. The differences are significant and maximal 
at ~ 100 km s" 1 , in the region of the spot. Figure 
[T3] shows the flux in these three emission lines as a 
function of distance from the QSO. It shows that, 
in the region where the radial velocities disagree, 
the flux distribution is very different in the three 
lines. Ha is maximal and [OIII] minimal at about 
1.5 from the QSO, i.e. close to the spot. [Oil] has 
intermediate strength. This is typical of a star for- 
mation region. Note that the absolute magnitude 
of the spot, as determined from the HST image, 
amounts to My ~ —17.3, which is als o rather typ- 
ical of starburst galaxies (Cox 19991) . The spot 
might thus correspond to the core of an intense 
star-forming region. 

Moving further away from the QSO, the ra- 
dial velocity curves tend to converge to roughly 
the same value, v ~ +120 km s _1 with respect to 
the QSO. The narrow emission lines become quite 
weak in this region which, in its outermost part 
(labelled F' in FigH2"]). also displays a broad (re- 
flected) Ha line. Combining the radial velocity 



shifts measured on the narrow lines and on the 
reflected broad line indicates that the reflecting 
medium is moving away from the QSO at a speed 
v ~ 370 km s -1 , and at an angle ~ 70° with 
respect to the line-of-sight. 

We are far from fully understanding the ori- 
gin of all the peculiarities found in this QSO host 
galaxy. The discrepant velocities measured from 
different narrow emission lines may be due to the 
fact that several gas clouds with different ioniza- 
tion stages contribute to the spectral emission in 
some regions of the host. In other regions, their 
velocities are sufficiently different for the different 
components to be resolved. Moreover, the strange 
behaviour of the radial velocity curve in the region 
of the spot, as well as the velocity shift between 
the QSO and the galactic bulge, indicate that this 
system is undergoing violent phenomena, which 
might be related to a collision between two galax- 
ies whose images are superimposed on the plane 
of the sky. 

6.5. HE 0354 5500 

The HST image (Fig. [5]) reveals a distorted 
host galaxy, with the probable presence of a few 
spiral arms, as well as an extended structure to the 
North, suggestive of a (nearly edge-on?) galaxy in 
violent interaction with the host. 

The diagnostic diagrams (Fig. [13)) show that 
the source of ionization changes as one moves 
along the slit, from the S-E (region A) to the N- 
W (region E) with a more important contribution 
from stellar light on the S-E parts (opposite to the 
colliding galaxy), a mixture of sources in the cen- 
tral regions and ionization completely dominated 
by the AGN in the N-W region E. 

Figure [3] shows that region E may contain some 
contribution from the two colliding galaxies. How- 
ever, its spectrum displays only emission lines, 
with no detectable continuum. Moreover, these 
emission lines extend to the outer regions, where 
no flux is measured on the HST image. We thus 
conclude that the spectrum is completely domi- 
nated by emission from AGN-ionized gas. The ra- 
dial velocity curve shows velocities increasing from 
+300 to +550 km s , as one moves away from the 
center. 

Moreover, the spectrum of this region E dis- 
plays a clear broad component in the Ha emis- 
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sion (Fig. [7]). This may be interpreted as ei- 
ther reflection of the QSO light by dust, or scat- 
tering by free electrons. The relative velocity of 
this broad component with respect to the QSO is 
v = +460 km s~\ which is close to the average 
radial velocity measured on the narrow emission 
lines. This near equality shows that the gas moves 
away from the central QSO, in a direction oppo- 
site to us, and at an angle ~ 45° with respect to 
the line-of-sight. 

Combining these two pieces of informations 
leads to the plausible interpretation that region 
E contains gas blown away by the AGN radiation. 
The radiation pressure, acting continuously as the 
gas moves away from the center, produces a con- 
tinuous acceleration and a speed increasing from 
about 400 to 600 km s" 1 in the region where we 
can measure it (i.e. outside the visible image of 
the host galaxy). 

This expanding cloud of ionized gas also scat- 
ters the QSO light. As it displays a range of ve- 
locities, the electron temperature, measured from 
the broadening of the Ha line (Table [2]), must be 
considered as an upper limit. However, this range 
of velocities (~ 200 km s _1 ) is too small to suffi- 
ciently broaden the line, as ~ 2000 km s _1 would 
be required. 

Other interesting features appear in the bulge 
region. The H/3 line displays a complex structure 
(Fig. [T5l top right), which can be interpreted as 
the sum of several components moving at different 
speeds. Indeed, the deconvolved HST image (Fig. 
IT5| left) shows, next to the bulge, three bright re- 
gions which, because of the spread of light due to 
the seeing, cannot be separated in the VLT spec- 
tra. Figure IT51 shows that the H/3 line can be fitted 
by 4 gaussians, which is also the number of emis- 
sion regions seen on the HST images. 

Due to the presence of the [Nil] lines, the Ha 
profile is not easy to interpret. However, Fig. [15] 
shows that a model with the same components as 
detected in H/3 can provide a satisfactory fit of the 
observations. 

On the other hand, only the strongest compo- 
nent is detected in the [OIII] lines. As this compo- 
nent is at the same radial velocity as the QSO, it 
probably corresponds to the bulge. Its position in 
the diagnostic diagrams (region C) shows a mix- 
ture of ionization sources, with a dominant contri- 



bution by shocks. On the other hand, the other 
components, with no detectable [OIII], should be 
dominated by stellar ionization. They probably 
correspond to starbursts initiated by the galactic 
collision, which might also be the cause of the QSO 
ignition. 

6.6. HE 0450 2958 

This very special case (Fi g. [T6"|) has been ana- 



lyzed in lMagain et al.l (|2005l ). No host is detected 
on the HST images, which allow to set an upper 
limit implying that it is underluminous by a factor 
of at least 6 fo r a QS O of that absolute magnitude. 
Merritt et all (|2006l) . noting that the QSO spec- 
trum displays rather narrow broad emission lines, 
suggest that the black hole might be of rather low 
mass, but with a high accretion rate. 

A companion galaxy is found at about 7 kpc 
S-E of the QSO. It has a strongly distorted shape, 
probably resulting from a collision. Careful re- 
moval of the QSO light reveals a compact emis- 
sion region (C) right next to the QSO. The size 
of this object is ~ lkpc and its spectrum shows 
that it consists of gas ionized by the QSO, with 
no detectable stellar continuum. 

The diagnostic diagrams in Fig. \W\ (bottom) 
show the clear separation between the companion 
galaxy spectrum and the other regions. The com- 
panion galaxy ionization is due to a mixture of 
star formation and shocks, whereas the other re- 
gions appear to be ionized by the QSO itself. This 
result tends to favour the view that, if there is star 
formation near the QSO, it is much weaker than 
the contribution of the ionization by the AGN. 

A broad Ha component is detected in the B 
region (Fig. [P7|) . in between the QSO and the 
companion galaxy. Its shift relative to the QSO 
Ha line center corresponds to a relative velocity 
of +745 km s^ 1 , the highest one measured in the 
present study. It has no relation cither with the ve- 
locities measured in the companion galaxy, which 
range from —60 to +200 km s -1 relative to the 
QSO. It is also completely different from the radial 
velocity measured from the narrow emission lines 
(+160 km s _1 ). Assuming that the matter respon- 
sible for the narrow line emissions and for the re- 
flection (or scattering) of the QSO light share the 
same motion (it may indeed belong to the same 
gas cloud), we can combine the two measured ra- 
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dial velocities to derive a speed of ~ +900 km s~ 
and an angle of 9 ~ 80°. Thus, this gas moves 
away from the QSO and towards the companion 
galaxy, in a direction close to the plane of the sky 
and at a speed close to 1000 km s _1 . This is one 
more indication that the system is undergoing vio- 
lent interactions and that this powerful QSO has a 
profound influence on its large-scale environment. 

7. General conclusions 

The present study has shown that: 

1. A joint analysis of spatially resolved VLT 
spectra and high-resolution HST images, 
both carefully processed with the MCS 
method, is a powerful way to characterize 
the QSO host galaxies; it allows to discover 
phenomena which would not be accessible 
through imaging or spectroscopy alone. 

2. The Ha line, in regions far away from the 
center, has a broad component due to re- 
flection of the nuclear light by dust or scat- 
tering by free electrons. Our study provides 
an extension of this phenomenon, only ob- 
served for Type 2 QSOs (i.e. objects in which 
the central engine is obscured by dust), to 
Type 1 QSOs. Evidence for extended scat- 
tered light has been independently found by 
Borguet et all (|2007t ). 



In most cases, the AGN has long-range ef- 
fects on its environment. This important 
point in the AGN-starbust interaction, al- 
ready known for bot h Type 1 (w ith the 



"proximity effects" , s eelBailtik at al.1 (119881) . 
Ferrarese fc Merrittl (E000h:ISiiaki et al.l(|2007l^ 



and Type 2 QSOs (jKauffman et alj|2003 ). is 
confirmed here in a completely different and 
independent way for Type 1 QSOs. 

We now summarize the main results on the in- 
dividual objects, in order of increasing AGN lumi- 
nosity. 

The two lower luminosity objects, HE 1503+0228 
and HE 2345-2906, which lie close to the some- 
what arbitrary limit between Seyfert galaxies and 
QSOs, are spirals. One of these spiral hosts has 
a very prominent bar and the other one has an 
extended arm or shell of gas and stars, probably 
reminiscent of a past interaction. 



HE 1434-1600 consists in a massive elliptical 
in gravitational interaction with a rather modest 
neighbor. Moreover, it contains shells of ionized 
gas blown up by the central QSO. 

HE 0306-3301, which has about the same lumi- 
nosity as the previous one, looks like a spiral but 
has an extremely odd radial velocity curve. The 
central QSO has a radial velocity of ~ 130km/s 
with respect to the galactic bulge. Moreover, a 
compact spot is found at ~ 6 kpc from the center. 
This spot, which might correspond to an intense 
starburst region, seems to add further perturba- 
tions to the radial velocity curve. 

HE 0354-5500 looks like a (barred ?) spiral 
undergoing a violent collision. It shows compact 
starburst regions close to the center and an ex- 
tended gas component most probably accelerated 
by the QSO radiation, up to distances as large as 
8 kpc. 

Finally, the most luminous object, HE 0450- 
2958 (which is a powerful infrared emittor and 
may also be the youngest QSO in the sample) 
has no detected host galaxy, but has a strongly 
distorted galactic companion and has long range 
influence on its neighborhood. 

Although our sample is quite small, the fact 
that the most powerful QSOs seem associated 
to the most dynamically perturbed host galaxies 
might not be a mere coincidence and we might 
be tempted to conclude that dynamical pertur- 
bations are efficient means of feeding the central 
AGN. Of course, even if the idea of a link between 
AGN activity and intera ctions is already 20 years 
old (panders et al.lll988l ). larger samples of bright 
radio-quiet QSOs with combined high resolution 
imaging and 3D spectroscopy are needed to fur- 
ther check these findings, to shed more light on 
the dynamical processes responsible for the igni- 
tion and fueling of the most active AGNs and to 
allow a more precise assessment of their impact on 
their galactic-scale environment. 

This work was supported by PRODEX Exper- 
iment Agreement 90195 (ESA and PPS Science 
Policy, Belgium). 
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Fig. 1.— Flatfielded observations of HE0354-5500 and a PSF star with HST/ACS. Top left: QSO long 
exposure. Top right: QSO short exposure. Bottom left: PSF star long exposure. Bottom right: PSF star 
short exposure. 
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Fig. 2. — Steps for QSO/host separation. 1) One of the long exposures of HE0354-5500 after reduction. 2) 
the point-like source has been removed from the observation. The PSF spikes are well removed. 3 to 6 are 
the results of the fit of the galaxy: 3) the analytical model found for the host galaxy bulge. 4) Numerical 
background. 5) The total fit of the galaxy (panel 3 + panel 4), the one used for analysis. Frame 6 displays 
the residuals (raw data minus total fit) of the whole process. 
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Fig. 3. Left: The VLT slit overlaid on the deconvolved image of the HE 0354-5500 host galaxy. Each 
column, corresponding to one pixel along the spatial direction of the slit, has a width of 0.2" and a height of 
1". The deconvolved HST image is used to identify the pieces of spectra which correspond to specific parts 
of the host galaxy. Right: The spatially resolved spectra are superimposed on the deconvolved host galaxy, 
which is here reconvolved by a Moffat in order to reach the VLT spatial resolution. 
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Fig. 4. — Spatially resolved spectra of several objects around H/3-[OIII]. Each column scopes different regions 
of the QSO environment. These regions are labeled A-E and defined in Fig. 1111141 and [TBI 
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Fig. 5. — Spectral modelling of the the D region in HE 2345 — 2906. The modeling of the spectrum by three 
gaussians and a linear continuum gives a very good fit. 




Fig. 6.— The fit of the QSO Ha line of HE 0354-5500 by two gaussians. 
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Fig. 7. — Left: 3-gaussian fit of the E region around the Ha line in HE 0354-5500, along with its residuals. 
Right: 5-gaussian fit, and residuals. Including a broad line with the QSO characteristics considerably 
improves the total fit of the Ha— [Nil] region. 
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Fig. 8. — HE 1503+0208. a) Reduced image, b) Deconvolved host galaxy with the different slits used to 
extract the spatially resolved spectra in overlay. The image reveals a regular spiral host surrounded by a 
kind of halo, which could be wider spiral arms, as suggested by the spectral analysis, c) Radial velocity 
curve, d) Diagnostic diagrams. 
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Fig. 9. — HE 1503+0228: The 5-gaussians fit of the outer region (top image) is improved when a stellar 
absorption component is added (bottom image). 
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Fig. 10. — HE2345-2906. a) Reduced image, b) Deconvolved host galaxy with the slits in overlay. The host 
is clearly a barred spiral galaxy. Wider arms can be found w 8kpc from the center, c) Radial velocity curve, 
d) Diagnostic diagrams. 
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Fig. 11. — HE1434-1600. a) Reduced image with the axe of the two radio lobes ([Condon et alJ Il998f ) 
superimposed, b) Deconvolved host galaxy with the slits in overlay, revealing filaments on both side of the 
center, and an undisturbed neighboring elliptical galaxy, c) The radial velocity curve tells us about the 
motion of the gaseous filaments, d) Diagnostic diagrams. 
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Fig. 12. — HE 0306-3301. a) Reduced image, b) Deconvolved host galaxy with the slits in overlay. As 
different smoothings had to be applied to the central regions and to the spot, a close-up of the spot with the 
appropriate smoothing is shown (with a point source), c) Radial velocity curve. In the region around — 1", 
the emission lines are resolved into 2 components. For Ha, the velocities of the 2 components are indicated. 
For the other emission lines, they are summarized by the crosses whose height corresponds to the range 
spanned by the measured velocities. The noise in the curve and the few discrepant points are due to the 
weakness of the emission lines, d) Diagnostic diagrams. The different ionizations measured on the separated 
components in the B-C region are indicated by open triangles. 
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Fig. 13. — HE 0306-3301: measured flux in the host galaxy spectrum for different emission lines ([OIII], 
[Oil] and Ha) as a function of distance from the QSO. The region E corresponding to the spot is indicated. 
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Fig. 14. — HE0354-5500. a) Reduced image, b) Deconvolved image with the slits in overlay (a zoom on the 
bulge is presented in Fig fT5|) . c) Radial velocity curve, d) Diagnostic diagrams. 
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Fig. 15. — Left: deconvolved central region of HE 0354-5500 revealing different bright areas. Right: fit of 
the spectrum of the central region C (see FiglM]) around H/3 (top panel) and Ha(bottom panel). 
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Fig. 16. — HE 0450-2958. a) Reduced image, b) Deconvolved image with the nucleus and the star removed 
thanks to the deconvolution process, c) Diagnostic diagrams. 
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Fig. 17. — Top image: fit of the region B in HE 0450-2958 around the Ha line. Bottom image: the different 
lines fit. The broad Ha line is strongly shifted compared to the narrow line. 
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